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A  Computer  Program  for 

Calculating  Inbreeding  Coefficients 

in  Dioecious  Populations 


By  Phillip  C.  Lowe,  research  poultry  husbandman,  Poultry  Research  Branch, 
Animal  Science  Research  Division,  Agricultural  Research  Service 


SUMMARY 

An  average  inbreeding  coefficient  for  each 
successive  generation  of  a  pedigree  dioecious 
population  is  obtained  by  calculating  a  weighted 
average  of  family  inbreeding  values.  The  source 
program  of  about  275  FORTRAN  IV  statements 
is  designed  for  punched  card  input  with  identi- 
fication of  sire,  dam,  and  progeny  on  each  record. 
One  generation  is  processed  with  each  computer 
run.  Computer  output  yields  data  on  covariances 
within  and  between  all  full-sib  families,  inbreeding 
coefficients  for  all  full-sib  families,  and  an  average 
inbreeding  value  for  the  current  generation. 

The  CDC  6500  computer  with  150,000  K 
(octal)  word  storage  used  for  program  development 
permitted  a  maximum  population  size  of  approxi- 
mately 200  dam  families  limited  to  no  more  than 
nine  progency  per  family.  The  number  of  families 
related  to  a  given  family  is  arbitrarily  limited  to 
99.  An  approximation  procedure  is  optional  if  this 
limit  is  exceeded. 

INTRODUCTION 

The  level  of  inbreeding  in  a  population  or  the 
rate  of  change  in  that  level  is  of  major  concern  to 
the  breeder,  geneticist,  or  researcher  since  progress 
in  breeding  programs  and  selection  systems  is 
usually  confounded  with  the  effects  of  inbreeding. 
Where  pedigree  records  of  a  population  are  main- 
tained, the  coefficients  of  inbreeding  of  individuals, 
families,  or  generations  may  be  calculated  by 
various  methods,  such  as  path  coefficients  or  co- 
variance  tables.  Hand  calculation  of  inbreeding 
coefficients  becomes  impractical  as  population  size 
and  number  of  generations  increase. 


Rehfeld  and  coworkers  (1967)  1  present  a  path- 
analysis  method  of  computing  inbreeding  by  a 
computer  program  adapted  to  a  moderate  memory 
capacity  computer  requiring  magnetic  tape  or  disk 
file  systems.  The  covariance  matrix  method  pre- 
sented herein  is  adapted  to  a  very  high  speed  com- 
puter with  large  central  memory  capacity.  All 
input,  output,  and  processing  is  accomplished  with 
punched  cards  and  central  memory;  this  eliminates 
the  need  for  magnetic  tape  or  disk  files. 

Wright  ( 1921)  2  showed  that  F=       , .       ,  where 

Fis  the  inbreeding  coefficient  of  the  individual  and 
COV  SD  the  covariance  between  the  individual's 
sire  and  dam.  This  permits  a  method  of  indexing 
all  possible  relationships  among  individuals  in  a 
given  generation  of  a  pedigreed  population  through 
a  table  or  matrix  of  such  covariances.  For  con- 
servation of  computer  running  time  and  memory 
storage,  the  covariances  are  best  calculated  and 
stored  on  a  family  basis  rather  than  on  an  indi- 
vidual basis. 

Only  three  types  of  covariances  need  be  used: 
(1)  covariances  between  members  of  different 
families,  (2)  covariances  between  members  of  the 
same  family,  and  (3)  covariances  of  any  individual 
of  a  family  with  itself  (this  equals  1  +  F,  where 
F  is  the  inbreeding  coefficient  for  each  member  of 
that  particular  family) . 

Since  all  covariances  between  members  of  a  full- 
sib  family  are  equal,  the  average  inbreeding  coeffi- 


1  Rehfeld,  C.  E.,  Bacus,  J.  W.,  Pagels,  J.  A.,  and  Dipert, 
M.  H.  Computer  calculation  of  Wright's  inbreeding  co- 
efficient, Jour.  Hered.  58:  81-84.  1967. 

2  Wright,  Sewall.  Systems  of  mating.  I-V.  Genetics  6: 
111-178.  1921. 


cient  for  a  given  generation  may  then  be  denned 
as  a  weighted  average  of  family  inbreeding  coeffi- 
cients calculated  as: 


i=  1 


where : 

-F\=the  inbreeding  coefficient  for  the  ith  family; 
7ij=the  number  of  individuals  in  the  ith  family 
(i—l  .  .  .  _/");  Fg=ihe  average  inbreeding  coeffi- 
cient for  a  generation. 

For  purposes  of  clarification,  an  example  of  a 
simple  five-generation  pedigree  with  punched  card 
input  and  output  and  printed  output  appears  on 
page  4.  Following  the  example  are  sections  con- 
taining the  flow  chart,  glossary  of  fortran  variables, 
and  the  source  program. 

DESCRIPTION  OF  PROGRAM 

The  source  program,  consisting  of  approxi- 
mately 275  FORTRAN  IV  statements,3  is  de- 
signed to  read  and  process  pedigree  information 
from  punched  cards.  Individual  animals  are 
assigned  family  codes  by  the  program  according 
to  their  full-sib  families.  Relationships  among  all 
coded  families  are  determined  and  included  in  the 
printed  and  punched  output.  As  subsequent 
generations  are  processed,  one  generation  per 
computer  run,  the  punched  output  from  the  pre- 
vious generation  containing  coded  families  and 
covariances  between  them  must  also  be  read  in 
and  stored.  The  relationship  between  any  two 
coded  families  of  the  current  generation  is  then 
determined  by  the  relationships  among  their 
parents. 

The  maximum  size  population  that  can  be 
processed  is  limited  by  the  computer's  storage 
capacity,  number  of  generations  to  be  processed, 
and  numbers  of  full-sib  families  in  a  given  genera- 
tion. Large  computers  with  storage  capacities  ot 
50  to  55K  words  accommodate  populations  having 
as  many  as  200  full-sib  families.  For  smaller  com- 
puters, the  arrays  specified  in  the  DIMENSION 
statement  and  listed  in  the  glossary  section  must 


3  Control  Data  6400/6500/6600  Computer  Systems, 
FORTRAN  Reference  Manual,  Publication  No. 
601749008.  1967. 


be  decreased  in  size  to  stay  within  the  limits  of  the 
computer's  storage  capacity.  Central  memory  and 
running  time  requirement  may  also  be  materially 
reduced  by  the  use  of  effective  parents  only  ( those 
parents  represented  by  their  progeny  in  the 
breeding  nucleus  of  the  next  generation) ;  this 
enables  larger  populations  to  be  processed  by  a 
computer  of  given  size.  The  program  arbitrarily 
limits  the  maximum  number  of  families  related  to 
a  given  family  to  99  and  the  maximum  number  of 
dams  mated  to  a  given  sire  to  nine. 

If  either  limit  is  exceeded,  an  error  message  is 
printed.  To  avoid  exceeding  the  maximum  num- 
ber of  related  families,  the  user  may  specify  a 
minimum  co variance  between  families  below  which 
families  are  assumed  to  be  unrelated.  An  average 
of  these  small  excluded  covariances  is  obtained  and 
appears  in  the  summary  of  the  printed  output.  A 
simple  method  is  provided,  as  explained  under 
"Preparation  of  parameter  cards,"  for  inputting  this 
value  into  the  subsequent  generation  where  it  is 
used  where  required,  as  part  of  the  inbreeding 
coefficient  of  families  in  that  generation.  All 
covariances  that  are  zero  or  are  smaller  than  the 
covariance  default  value  specified  in  parameter 
card  No.  1  are  omitted  from  the  covariance  matrix. 
This  procedure,  which  is  designed  to  reduce  com- 
puter storage  requirements,  permits  the  use  of  two 
rectangular  two  dimensional  arrays,  IR  and  IC, 
for  storage  of  the  identification  of  related  families 
and  their  covariance  values,  respectively. 

PREPARATION  OF 
PARAMETER  CARDS 

1.  Parameter  card  No.  1. 

Columns  Description 

1-3       Generation  number 
4-6       Number  of  cards  in  pedigree  deck 
7-9       Estimated  maximum  number  of  families 
related  to  a  given  family — limit  set  at  99 
in  program 
10-12     Estimated  maximum   number  of  full-sib 

families 
13-15     Estimated   maximum   number   of   full-sib 

families  plus  1 
16-19     Covariance   default   level — minimum   de- 
sired covariance  multiplied  by  4 
20-23     Average  value  of  excluded  covariances  in 
previous  generation 


2.  Parameter  card  No.  2.  Columns  1  to  6  equal 
IBLANK,  a  data  word  used  to  check  on  the 
number  of  cards  contained  in  the  pedigree  deck. 

3.  Parameter  card  No.  3  is  a  blank  trailer  card. 

PREPARATION  OF  PEDIGREE 
AND  DATA  DECKS 

1.  A  pedigree  deck  containing  identification  of 
sire,  dam,  and  progeny  in  three  four-digit 
identification  fields  is  prepared  for  each  genera- 
tion. One  punched  card  is  required  for  each 
individual  progeny  in  the  population.  The 
source  program,  as  listed  on  p.  14,  is  written 
for  a  data  format  that  provides  identification 
for  sire  in  columns  5  to  8  and  identification  for 
dam  in  columns  11  to  14  and  identification 
for  progeny  in  columns  17  to  20.  Alternative 
data  formats  in  the  pedigree  deck  may  easily 
be  accommodated  by  appropriate  changes  in 
FORMAT  statement  No.  104.  The  pedigree  deck 
must  be  sorted  by  dams  within  sires  in  ascend- 
ing order.  Comment  cards  and  explanatory 
statements  to  assist  the  user  in  preparation  of 
input  decks  appear  at  frequent  intervals  in  the 
source  program  and  in  the  printed  computer 
output. 

2.  A  data  deck  containing  family  codes,  a 
covariance  matrix,  and  family  inbreeding 
coefficients  is  produced  by  the  computer  as  a 
part  of  the  output  after  the  processing  of  each 
generation.  This  deck  is  regarded  as  output 
relative  to  the  generation  currently  being 
processed.  But  it  must  be  rearranged  and 
prepared  as  input  to  accompany  the  new 
pedigree  deck  during  processing  of  the  subse- 
quent generation. 

The  punched  output  (data  deck)  for  each  gen- 
eration consists  of  three  sections  each  separated 
by  a  blank  card.  To  rearrange  this  output  as 
input  for  the  next  generation:  (1)  Remove  the 
last  section  of  usually  25  cards  or  less  and  discard 
the  blank  card  preceding  it;  (2)  remove  the 
blank  card  separating  section  one  (30  cards  or 
less)  and  section  two  and  insert  section  three  at 
that  location  in  the  deck.  No  data  deck  is  neces- 
sary for  the  processing  of  the  first  generation. 
The  order  of  input  decks  is  outlined  below. 
First  generation 

1.  Source  deck 

2.  Parameter  card  No.  1 

3.  Pedigree  deck 


4.  Parameter  card  No.  2 
Second  and  later  generations 

1.  Source  deck 

2.  Parameter  card  No.  1 

3.  Pedigree  deck 

4.  Parameter  card  No.  2 

5.  Data  deck — reorganized  as  discussed  above 

6.  Parameter  card  No.  3 


ERROR  MESSAGES 

The  program  provides  for  four  types  of  error 
messages  as  follows: 

1.  PARAMETER  CARD  NO.  2  OMITTED  OR 
INCORRECT  VALUE  FOR  NUMBER  OF 
CARDS  IN  PEDIGREE  DECK 

This  message  is  printed  during  processing  time 
if  parameter  card  No.  2  is  erroneously  placed 
or  omitted  or  when  the  incorrect  value  for  the 
number  of  cards  in  the  pedigree  deck  is  punched 
into  columns  4  to  6  of  parameter  card  No.  1. 
Processing  is  terminated  when  this  condition 
is  encountered.  To  correct  this  condition, 
recount  cards  in  pedigree  deck  and  check  for 
presence  of  parameter  card  No.  2. 

2.  INVALID  PEDIGREE— PREVIOUS  GEN- 
ERATION PROGENY  CODE  OF  S1RE= 
ss,  OF  BAM=dd,  CURRENT  GENERA- 
TION PROGENY  CODE  OF  INDIVID- 
UAL =ii 

This  message  is  printed  during  processing  time 
if  a  search  of  the  covariance  matrix  reveals  that 
either  or  both  parents  of  an  individual  in  the 
current  generation  were  not  present  in  the 
previous  generation.  (The  notations  ss  and  dd 
represent  the  sire's  and  dam's  family  code, 
respectively,  if  they  are  present  in  the  previous 
generation,  otherwise  the  value  will  be  zero. 
The  symbol  ii  represents  the  progeny  code  of  an 
individual  in  the  current  generation.)  After  this 
interruption,  processing  continues  with  the 
next  family.  The  suggested  remed}7  is  to  correct 
the  pedigree  error  and  rerun  the  previous 
generation. 

3.  NUMBER  OF  RELATED  FAMILIES  EX- 
CEEDS xx,  PROGENY  CODE  OF  NEXT 
RELATED  FAMILY =cc 

The  message  occurs  when  the  number  of 
families  related  to  a  given  family  exceeds  the 
value  placed  in  columns  7  to  9  of  parameter 


card  No.  1.  Here  xx  equals  the  maximum 
number  of  related  families,  and  cc  is  the 
progeny  code  of  the  next  family  in  the  covar- 
iance  matrix.  The  problem  is  eliminated  by 
increasing  the  magnitude  of  the  covariance 
default  value.  Processing  is  not  interrupted. 
SIRE  xx  MATED  TO  MORE  THAN  9 
DAMS 

The  statement  occurs  when  there  are  more  than 
nine  full-sib  families  within  one  sire  family. 
The  sire  is  identified  by  the  value  xx,  which 
is  the  sequential  number  of  the  sire  in  the 
pedigree  of  the  current  generation.  Processing 
is  terminated. 


EXAMPLE  OF  CARD  INPUT  AND 

OUTPUT  FOR  A 
FIVE-GENERATION  POPULATION 

Identification  Nos.  for  a  five-generation  pedigree 

population  used  in  the  example  on  the  following 

pages  are  as  follows: 

Identification  Nos. 


Input  for  First  Generation 


Generation 
1 


Sire 
115 
117 
163 


3d 


1531 


1535 

1540 
1626 


Dam 
2528 
2617 
8020 


23 


24 

1536 


1541 

1622 
1542 
1543 


Progeny 

163 

8020 

23 

24 

30 

1531 

1535 

1536 

1541 

1540 

1542 

1543 

1622 

1626 

23 

15 

12 


Program 

00100200200200300000000 
115  2528  163 
117  2617  8020 

IBLANK 

EOF 


(Source  Deck) 
(Param.  Card  No.  1) 
(Pedigree  Deck) 

(Param.  Card  No.  2) 
(End  of  File  Card) 


Punch  Output  From  First  Generation 

2    PROGENY   CODES   DERIVED    FROM    2 
FAMILIES  OF  FIRST  GEN. 

163  118020  21 
Blank  Card 
1000010000 
1  11  5000 

1  21  5000 

Input  for  Second  Generation 

Program 

00200300200200300000000 
163     8020       23 

24 

30 
IBLANK 

2  PROGENY  CODES  DERD7ED 
FROM  2  FAMILIES  OF  FIRST 

GEN.  I  DATA 

163     118020     21  I      DECK 

1000010000 
1   11  5000 
1  21   5000 
Blank  Card     (Param.  Card  No.  3) 
EOF 

Punch  Output  From  Second 
Generation 

3  PROGENY  CODES  DERIVED  FROM 
1  FAMILY  OF  SECOND  GEN. 
23     11     24     11     30     11 
Blank  card 
1  11  5000 
Blank  card 
10000 


Input  for  Third  Generation 

Program 

00300400200200300000000 

30    23  1531 

1535 

1536 

24  1541 

IBLANK 

3  PROGENY  CODES  DERIVED  FROM 

1  FAMILY  OF  SECOND  GEN. 
23     11     24     11     30     11 

10000 

1  11     5000 
Blank  card 
EOF 

Punch  Output  From  Third  Generation 

4  PROGENY  CODES  DERIVED  FROM 

2  FAMILIES  OF  THIRD  GEN. 
1531     111535     111536     111541     12 
Blank  card 

2  11  7500  12  6250 
2  12  7500  11  6250 

Blank  card 
12500  12500 

Input  for  Fourth  Generation 

Program 

00400500200200300000000 
1531     1536     1540 


1535  1541 


1542 
1543 
1622 
1626 


IBLANK 

4  PROGENY  CODES  DERIVED  FROM 
2  FAMILIES  OF  THIRD  GEN. 
1531     111535     111536     111541     12 
1250012500 
2   11  7500  12  6250 
2   12  7500  11  6250 
Blank  card 
EOF 


Punch  Output  From  Fourth 
Generation 

5    PROGENY    CODES   DERIVED    FROM    2 
FAMILIES  OF  FOURTH  GEN. 

1540     111542     111543     111622     211626     21 
Blank  card 

2  1110000  21  6875 

2  21  9375  11   6875 
Blank  card 
1375013125 


Input  for  Fifth  Generation 


Program 
00500300300300400000000 


1540  1622 

1626  1542 

1543 

IBLANK 

5    PROGENY 


23 

15 
12 


FROM    2 


CODES   DERIVED 
FAMILIES  OF  FOURTH  GEN. 

1540  111542  111543  111622  211626  21 
1375013125 

2  1110000  21  6875 

2  21  9375  11  6875 
Blank  card 
EOF 


Punch  Output  From  Fifth  Generation 

3    PROGENY   CODES   DERIVED    FROM    3 
FAMILIES  OF  FIFTH  GEN. 
23     11      15     21      12     22 
Blank  card 
3  1110156  21  8281  22  8281 
3  2110156  11  8281  22  9218 
3  2210156  11  8281   21  9218 
Blank  card 
134371343713437 


FACSIMILE  OF  COMPUTER 

PRINTOUT  OF 

FIVE-GENERATION  EXAMPLE 

Printed  Output  From  First  Generation 

2    PROGENY   CODES   DERIVED    FROM. 2 

FAMILIES  OF  FIRST  GEN. 
CODE     CODE     CODE     CODE    CODE 
11      163     21    8020 

CURRENT  GENERATION  RELATED  FAM- 
ILIES   AND    COVARIANCES    WITHIN 
AND  BETWEEN  FAMILIES  FOLLOW. 
NUMBER 
OF  REL.  FAM.     FAM.  COV.     FAM.  COV. 
1  11  5000 

1  21  5000 

THE  FOLLOWING  ARE  FAMILY  INBREED- 
ING COEFFICIENTS  FOR  THE  CUR- 
RENT GENERATION. 
.0000  .0000 


Printed  Output  From  Second 
Generation 

3  PROGENY  CODES  DERIVED  FROM  1 
FAMILY  OF  SECOND  GEN. 

CODE     CODE     CODE    CODE      CODE 

11        23      11        24     11       30 

CURRENT  GENERATION  =2 

NUMBER  OF  MASTER  CARDS=3 

PUNCH  OUTPUT  FROM  2  FAMILIES  IN 
PREVIOUS  GEN.  HAS  BEEN  READ  IN 
AND  STORED. 

CURRENT  GENERATION  RELATED  FAM- 
ILIES AND  COVARIANCES  WITHIN 
AND    BETWEEN   FAMILIES   FOLLOW. 

NUMBER  OF  REL.  FAM.  FAM.  COV.  FAM.  COV. 
1  11  5000 

THE  FOLLOWING  ARE  FAMILY  INBREED- 
ING COEFFICIENTS  FOR  THE  CUR- 
"  RENT  GENERATION. 

.0000 

.. -SUMMARY  FILE. ..SUMMARY  FILE___ 

GENERATION  NO.  =2 

NUMBER  OF  INDIVIDUALS  =3 

NUMBER  OF  SIRES  =1 

NUMBER  OF  DAMS  =1 

NUMBER  OF  INVALID  PED.=0 


AVERAGE  NUMBER  OF  RELATED  FAMI- 
LIES IN  THIS  GEN.  =1.0 

NUMBER  OF  DROPPED  COVARIANCES— 
CURRENT  GENERATION=0.0 

AVERAGE  DROPPED  COVARIANCE— PRE- 
VIOUS GENERATION         =0.0000 

AVERAGE  DROPPED  COVARIANCE— CUR- 
RENT GENERATION  =0.0000 

AVERAGE  PERCENT  INBREEDING— CUR- 
RENT GENERATION  =0.0000 

Printed  Output  From  Third 
Generation 

4    PROGENY   CODES   DERIVED   FROM    2 

FAMILIES  OF  THIRD  GEN. 
CODE    CODE    CODE    CODE 
11       1531     11       1535     11       1536     12       1541 

CURRENT  GENERATION  =3 

NUMBER  OF  MASTER  CARDS =4 

PUNCH  OUTPUT  FROM  1  FAMILY  IN 
PREVIOUS  GEN.  HAS  BEEN  READ  IN 
AND  STORED. 

CURRENT  GENERATION  RELATED  FAM- 
ILIES   AND    COVARIANCES    WITHIN 
AND   BETWEEN   FAMILIES   FOLLOW. 
NUMBER 

OF  REL.  FAM.   FAM.  COV.   FAM.  COV. 
2  11  7500       12  6250 

2  12  7500       11  6250 

THE  FOLLOWING  ARE  FAMILY  INBREED- 
ING COEFFICIENTS  FOR  THE  CURRENT 
GENERATION. 

.2500  .2500 

__ -SUMMARY  FILE  .-.SUMMARY  FILE_._ 

GENERATION  NO.  =3 

NUMBER  OF  INDIVIDUALS  =4 

NUMBER  OF  SIRES  =1 

NUMBER  OF  DAMS  =2 

NUMBER  OF  INVALID  PED.=0 

AVERAGE  NUMBER  OF  RELATED  FAMI- 
LIES IN  THIS  GEN.  =2.0 

NUMBER  OF  DROPPED  COVARIANCES- 
CURRENT  GENERATION=0.0 

AVERAGE  DROPPED  COVARIANCE— PRE- 
VIOUS GENERATION  =0.0000 

AVERAGE    DROPPED     COVARIANCE 
CURRENT  GENERATION  =  0.  0000 

AVERAGE  PERCENT  INBREEDING- 
CURRENT  GENERATION =25.  0000 


Printed  Output  From  Fourth 
Generation 

5    PROGENY    CODES   DERIVED    FROM   2 
FAMILIES  OF  FOURTH  GEN. 

CODE    CODE    CODE    CODE    CODE 

11    1540     11    1542     11   1543     21    1622     21   1626 

CURRENT  GENERATION  =4 

NUMBER  OF  MASTER  CARDS  =5 

PUNCH  OUTPUT  FROM  2  FAMILIES  IN 
PREVIOUS  GEN.  HAS  BEEN  READ  IN 
AND  STORED. 

CURRENT  GENERATION  RELATED  FAM- 
ILIES   AND    COVARIANCES    WITHIN 
AND   BETWEEN   FAMILIES   FOLLOW. 
NUMBER  OF 
REL.  FAM.       FAM.  COV.     FAM.  COV. 
2  1110000  21  6875 

2  21  9375  11  6875 

THE  FOLLOWING  ARE  FAMILY  IN- 
BREEDING COEFFICIENTS  FOR  THE 
CURRENT  GENERATION. 

.3750  .3125 

. ..SUMMARY  FILE.. .SUMMARY  FILE___ 

GENERATION  NO.  =4 

NUMBER  OF  INDIVIDUALS  =5 

NUMBER  OF  SIRES  =2 

NUMBER  OF  DAMS  =2 

NUMBER  OF  INVALID  PED.  =  0 

AVERAGE  NUMBER  OF  RE- 
LATED FAMILIES  IN  THIS 
GENERATION.  =2. 0 

NUMBER  OF  DROPPED  COVARI- 
ANCES—CURRENT  GENERA- 
TION. =  0.  0 

AVERAGE  DROPPED  COVARI- 
ANCE— PREVIOUS  GENERA- 
TION. =  0.  0000 

AVERAGE  DROPPED  COVARI- 
ANCE— CURRENT  GENERA- 
TION. =  0.  0000 

AVERAGE  PERCENT  INBREED- 
ING—CURRENT GENERA- 
TION. =  35.  0000 


Printed  Output  From  Fifth 
Generation 

3    PROGENY   CODES   DERIVED    FROM    3 
FAMILIES  OF  FIFTH  GEN. 
CODE    CODE    CODE 
11  23     21  15     22  12 

CODE    CODE    CODE 
CURRENT  GENERATION  =5 

NUMBER  OF  MASTER  CARDS=3 
PUNCH  OUTPUT  FROM   2   FAMILIES  IN 
PREVIOUS  GEN.  HAS  BEEN  READ  IN 
AND  STORED. 
CURRENT  GENERATION  RELATED  FAM- 
ILIES   AND    COVARIANCES    WITHIN 
AND   BETWEEN   FAMILIES   FOLLOW. 
NUMBER  OF 
REL. 

FAM.    FAM.  COV.  FAM.  COV.  FAM.  COV. 
3  1110156  21  8285  22  8281 

3  2110156  11  8281  22  9218 

3  2210156  11  8281  21  9218 

THE  FOLLOWING  ARE  FAMILY  INBREED- 
ING  COEFFICIENTS  FOR   THE   CUR- 
RENT GENERATION. 
.3437  .3437  .3437 

...SUMMARY  FILE. ..SUMMARY  FILE... 
GENERATION  NO  =5 

NUMBER  OF  INDIVIDUALS  =3 
NUMBER  OF  SIRES  =2 

NUMBER  OF  DAMS  =3 

NUMBER  OF  INVALID  PED.  =  0 
AVERAGE  NUMBER  OF  RELATED 

FAMILIES  IN  THIS  GEN.         =3.  0 
NUMBER  OF  DROPPED  COVARI- 
ANCES—CURRENT  GENERA- 
TION =  0.  0 
AVERAGE     DROPPED     COVARI- 
ANCE— PREVIOUS    GENERA- 
TION =  0.  0000 
AVERAGE     DROPPED     CO  VARI- 
ANCE—CURRENT    GENERA- 
TION =  0.  0000 
AVERAGE    PERCENT    INBREED- 
ING—CURRENT       GENERA- 
TION =  34.  3700 


GLOSSARY  OF  FORTRAN  VARIABLES 

Variable  Type  '  Definition 

Number  of  covariances  below  covariance  default  value 

Average  defaulted  covariance — previous  generation 

Average  defaulted  covariance — current  generation 

Average  percent  inbreeding — current  generation 

Average   number    of   defaulted    covariances    per   family   in    the    current 

generation 
Average  number  of    families  related    to  a    given  family  of    the  current 

generation 
Family  inbreeding  coefficients — current  generation 
Sum  of  covariances  among  the  parents  of  related  families 
Sum  of  defaulted  covariances — current  generation 
Covariance  default  value  (covariance  X  4) 
Covariance  between  any  two  families 
Number  of  individuals  over  which  inbreeding  is  averaged 
Number  of  families  in  the  current  generation 

Number  of  related  families  with  covariances  exceeding  the  default  value 
Stored  family  codes  of  previous  generation 
Progeny  codes  representing  families  of  the  current  generation 
DO  loop  index 

Value  defined  in  the  DATA  statement 
Dummy  variable 

Covariance  between  families  in  the  previous  generation 
Covariance  between  any  two  families  of  the  current  generation 
Identification  of  dam 
Part  of  family  code  designating  she  line 
Code  of  dam  family 
Code  of  sire  family 

Generation  number  of  current  generation 
DO  loop  index — number  of  related  families 
Covariances  within  families  of  the  current  generation 
Covariances  within  families  of  previous  generation 
Summation  of  weighted  family  inbreeding  values 
Identification  of  individual  or  progeny 
Number  of  progeny  per  full-sib  family 
Subscript  used  in  a  DO  loop 

An  array  of  families  related  to  a  given  family  of  the  previous  generation 
An  array  of  families  related  to  a  given  family  of  the  current  generation 
Identification  of  sire 
Individual's  identification   associated   with  family  code  of  the  previous 

generation 
Individual's   identification    associated    with    family   code   of   the   current 

generation 
DO  loop  subscript — progresses  from  1  to  NUMFPI 
DO  loop  index 

Variable  used  in  establishing  values  for  family  codes 
DO  loop  index — number  of  full-sib  families  in  current  generation 
Subscript  value 

See  footnote  at  end  of  glossary. 


ACCNUM 

R 

AVDROP 

R 

AVGDRP 

R 

AVGINB 

H 

AVNODR 

R 

AVNORL 

R 

COEFF 

R  (D,E) 

COV 

R 

COVACC 

R 

COVDEF 

R 

CV 

R 

D 

R 

Dl 

R 

D2 

R 

FCODEI 

I  (D) 

FCODEP 

1(D) 

I 

I 

IBLANK 

I  (DATA) 

IBLK 

I 

IC 

1(D) 

ICOV 

1(D) 

IDAM 

1(D) 

IDIO 

I 

IFAMD 

I  (D) 

IFAMS 

I  (D) 

I  GEN 

I 

INBR 

I 

INBRDC 

I  (D,  E) 

INBRDP 

1(D) 

INBSUM 

I 

IND 

1(D) 

INDWIN 

I  (D) 

IP1 

I 

IR 

I  (D) 

IREL 

I  (D) 

ISIRE 

1(D) 

ISTORI 

I  (D,  E) 

ISTORP 

I  (D,  E) 

11 

1 

J 

I 

JD10 

1 

JP 

I 

JPl 

I 

» 


Variable 


Type 


K 

I 

KACCUM 

I 

KCOUNT 

I 

L 

1(D) 

LIMIT 

I 

LIMITD 

I 

LIMITS 

I 

LISTED 

I 

LISTRS 

I 

LS 

I 

LT 

I 

LU 

I 

LVALUE 

I 

LI 

I 

M 

I 

MALES 

I 

MATCHD 

I 

MATCHS 

I 

N 

I 

ND 

I 

NEWFAM 

1(D) 

NFC 

I 

NFLD10 

I 

NFNDIO 

I 

NFP 

I 

NI 

I 

NIOLD 

I 

NIP 

I 

NIPl 

I 

NM 

I 

NMP1 

I 

NUMFAM 

I 

NUMFPl 

I 

NUMINV 

I 

NUMREL 

I 

TOTINB 

R 

X 

R 

Definition 

Counter — number  of  related  families 

Cumulative  number  of  related  families  per  generation 

DO  loop  parameter — number  of  full-sib  families  in  current  generation 

A  control  variable  assigned  to  number  of  related  families 

DO  loop  parameter — number  of  related  families 

DO  loop  parameter — number  of  families  related  to  female  parent's  family 

DO  loop  parameter — number  of  families  related  to  male  parent's  family 

DO  loop  index — number  of  co variances  among  families  related  to  dam 

DO  loop  index — number  of  covariances  among  families  related  to  sire 

DO  loop  index — number  of  full-sib  families  in  current  generation 

DO  loop  index — number  of  full-sib  families  in  previous  generation 

DO  loop  index — number  of  full-sib  families  in  previous  generation 

Implied  DO  loop  parameter — number  related  families 

Implied  DO  loop  index — number  related  families 

Implied  DO  loop  index — number  related  families 

Number  of  full-sib  families  in  current  generation 

Subscript — number  families  related  to  dam 

Subscript — number  families  related  to  sire 

DO  loop  index 

Number  of  families  in  covariance  table 

Family  currently  being  checked  through  covariance  table 

Number  of  full-sib  families  in  current  generation 

Test  variable  used  in  detecting  full-  or  half-sib  families 

Test  variable  used  in  detecting  full-  or  half-sib  families 

Number  full-sib  families  in  previous  generation 

Number  of  individuals  in  current  generation 

Number  of  individuals  in  previous  full-sib  family 

Number  of  individuals  in  previous  generation 

NIP+1 

Number  of  cards  in  pedigree  deck 

NM  +  1 

Estimated  maximum  number  of  full-sib  families  in  current  generation 

NUMFAM +  1 

Number  of  invalid  pedigrees 

Maximum  number  of  related  families  expected 

Sum  of  inbreeding  coefficients  in  all  full-sib  families 

Inbreeding  in  parents  of  a  given  family 


1 1,  Integer;  R,  real;  D,  dimension;  and  E,  equivalence. 
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E> 


PRINT 

COVARIANCES 

MATRIX  AND 

FAM.  INBR. 


READ  AND  STORE 

PROG.  CODES 

FROM 

PREVIOUS 

GEN. 


READ  AND  STORE 

INBR.  COEFF. 

FROM 

PREVIOUS 

GEN. 


ASSIGN  CODES 
TO  PARENTS  OF 

FAM  IN 
CURRENT  GEN. 


READ  AND  STORE 

COV.  MATRIX 

OF  PREV. 

GEN. 


F  (FAMILY)  = 
AVDROP 


NO 
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MATRIX  FOR 
PARENTS  OF 
EACH  FAM. 


F   (FAMILY)  ■■ 
1+    COVSD 

' 2 


NO 


INVALID 
PEDIGREE 


V 
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COMPUTE 

COVARIANCES 

WITHIN   FAMILIES 


/PARENTsX 
/        ARE        \— 

COV.  = 
(INBR. 

(COV.  + 
OF  SIRE 

NFULL-SIBS/' 

FAM.  +   INBR.  OF 
DAM   FAM.)  /2)  /2 

J  YES 

COV.  =  (COV.  + 
INBR.  OF  PARENT 

FAM.)  /2 

t 

COMPUTE 

COV.  BETWEEN 

FAMILIES 

SEARCH  COV. 

MATRIX  FOR 

RELATED  PARENTS 


NO 


PARTIAL  COV. 
COV.  +  COV. 


YES 


HALF-SIB  \       YES 
.RELATIONSHIP/ 


\? 


DAM 

^/related  tq\      yes      » 


PARTIAL  COV. 
COV.  +  INBR. 


PARTIAL  COV. 
COV.  +  COV. 
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PUNCH 
BLANK 
CARD 


PUNCH  AND  PRINT 

FAMILY 

INBREEDING 

VALUES 


PARTIAL  COV. 
COV.  +  COV. 


COV.  BETWEEN 
TWO  FAMILIES  = 

SUM  OF 
PARTIAL  COVS. 


YES 


NO 


ACCUMULATE 
DEFAULT 

YES  /  C0V     \ 

fES 

\  VALUE  / 

Yno 
/related\    . 

COVARIANCES 

PRINT  ERROR 
MESSAGE 

\     LIMIT     / 
f  NO 

COV.  =   (SUM  OF 
PARTIAL  COV.)  /4 

& 


PRINT  AND  PUNCH 

COV.  TABLE 

FOR  ONE 

FAMILY 


^ 
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PROGRAM  MAIN  ( INPUT i OUTPUT , PUNCH « 7 APFS  = INPUT ♦ TAPE6=0UTPUT« TAPF7sPU 
INCH) 
C 

C  PROGRAM  TO  COMPUTE!  AVG.  INRR.  COEFFICIENT  USING  FAMILY  PEMJGRFF*** 
C  PLACE  PARAME7FR  CARD  NO.  1  AHEAD  OF  PFDIGPFF  DECK««**»*»*«*»*»*»»* 
C  OUTPUT  FROM  PREVIOUS  GEK-.  IS  PART  nF  TNPUT  FO«  CURRENT  GEN. ******* 
C  PLACE  ****  IHLANK  ****CAKD  REHIND  LAST  CARD  OF  PEnlGREE  DECK****** 
C  PLACE  BLANK  CARD  BEHIND  LAST  CARD  OF  DATA  DECK******************** 
C  SORT  PEDIGREF  DECK  BY  DAM/SIRE  IN  ASCFNDTNG  ORDER**********'******* 
C 

EQUIVALENCE  < TSTuRP, ISTORj j  ,  (COEFF ♦ INRRDC) 
INTEGFR  FCODEPtFCODEI 

DIMENSION  FCODEPOOD*  lSTORP(30l>,  FCOREI(301),  TR(?05.9W)»  IC(?0 
15,99).  ISTORI(301),  IFAMS(?50),  IFAMr>(?50)t  NEWFAm(250),  ISIRE(301 
2)»  IDAMOOD*  INU(300)»  iRELtOR)*  TC0V(99)»  L(251)»  INDWIN(2B]),  T 
3NBRDP(251),  lNBRUC(251),  COEFF<2^1) 
DATA  IBLANK/6HIBLANK/ 
C 

C  READ  PARAME7EP  CARD  NO.  1  WITH  GEN.  NO., NO.  OF  PFnlGRFE  CARDS, ETC. 
C 

READ  (5»101)  IGEN,NM»NUHRFL»NHMEAM,N|JMFP1,C0VDEF,AVDR0P 
NUMINVsO 
KACCUM=0 
INBSUmsO 
COVaCc=0. 
ACCNUm=0. 
DO  10?  I=V,NUMREL 
IRFL(I)=0 
10?  ICOV(T)=0 
NMPlsMM+1 
DO  10^  I=I,NMPl 
ISTRE(I)=0 
IDAM(I)=0 
103  IND(I)bO 
C 

C      READ  TN  PFDIGREE  DECK  CONTAINING  PFDIORFE  IDENTIFICATION********** 
C      FORMAT  STATEMENT  NO  104  MAY  BF  CHANGED  TO  ACCOMODATE  VARIOUS  DATA 
C      FORMATS.* ************************************** ******************* 
C 

DO  10S  I=1,NM 

READ  (5*104)  ISIRE(I) «IDAM(T) »IND(T) 
10*;  CONTINUE 

READ  (5tlf)6)  TBLK 
IF  (IRLK.EO.IRLANK)  GO  TO  108 
WRI7E  (6tl07) 
GO  TO  191 
10R  DO  109  I=1,NUMFP1 

INRRDC(I)=10000 
10P    L(I)=6 

c 

C  GENERATE  FAMILY  CODES  AND  OUTPUT  FnR  USE  IN  NEXT  GEnFPAT ION******* 
C 

NFC  =  0 
NI  =  0 

DO  US  1  =  1,60 
DO  112  J=l,9 
110        Nl=Nl*l 

NIP1=NI*1 

FCODEP(NI)=I*10»J 

ISTORP(Nl)slND(NI) 


1  1 


IF  (IDAM(NI)-IDAM(NIPI))  111,110*111 

11]        NFCrNFC^l 

IF  (ISIRE(NI)-ISIRE(MIPI)  J  H4,H?,iU 
11?     CONTINUE 

WRITE  (6.H3)  ISIHE(NI) 

GO  TO  191 

114  IF  (ISTRE1NI*!))  115.116,115 

115  CONTINUE 

lift  PUNCH  117,  NI,NFC»IGEN 

PUNCH  118,  (ISTOHP(I) ,FCODEP(I) ,1=1 .NT) 
PUNCH  119 

WRITE  (6.120)  NI.NFC.IGEN 
WRITE  (6.121) 

WRITE  (6.122)  (FCODEP(I) .ISTORP(I)  .I  =  1,NI) 
MALES=FCODEP(NI)/10 
IF  (Ir,EN-l)  123.123.136 
C 

C      COMPUTE  COVAPIANCES  WITHIN  ANT  RETwEEN  FAMILIES  In  FIRST  (JFNER.»»» 
C 

123  WRITE  (6.124) 
WRITE  (6.152) 
DO  13?  1=1, NI 
IPl«!*l 
K  =  l 

IF  (FCOOEP(I) .EQ.FCOOEP(IPl) )  Go  TO  132 
DO  130  J=1.NI~ 
JP1=J*1 

IF  (FCODEP(J) .EQ.FCOOEP(JPI) )  GO  TO  130 
IF  (FCODEP(I)-FCOOEP(J) )  1?8,127,12R 
127        ICOV(1)=5000 

IREL(D=FCul)EP(I) 
GO  TO  130 
12B        ID10=FCODEP(I)/10 
JD10=FCODEP(J)/10 
IF  (ID10-JU10)  131.129,130 

129  KsK*l 
ICOV(K)=25U0 
IREI.(K)sFCUOEP(J) 

130  CONTINUE 

131  PUNCH  134.  K, (IREL(Ll) ,TC0V(L1> ,Ll=l,K) 
WRTTE  (6.126)  K»  (IR£I.  (M)  ,  ICOV  (M)  ,M=1  ,K) 

13?  CONTINUE 
PUNCH  119 

PUNCH  139,  (INBRUC(J) »J=1,nFC) 
DO  131  JPsl.NFC 

C0FFF(JP)=TNBH0C(JP) 
COFFF(jP)=(COtFF(JP)/lO000,)-l. 
137  CONTINUE 

WRITE  (6.187)  (COEFF(J) ,J=i,NFC) 
GO  TO  191 
C 

C      READ  TN  AND  STORt  FAMll*  COOES  FROM  PREVIOUS  GEN. »«»«»*»«•#»»•»#«* 
C 

13ft  READ  (5.137)  NIP.NFP 

READ  (5.138)  (ISTORI (I) ,FCOOEI (I) »T=1,NIP) 
C 

C      READ  TN  AND  STORt  INBREEDTNG  COEFFICIENTS  FROM  PREVIOUS  GENERATION 
C 

READ  (5.139)  (INbRDP(J) .J=1,NFP) 
C 
C      ASSIGN  CODES  TO  PARENTS  OF  FAMILIES  IN  CURRFNT  GEM.  SO  THAT  PARFMT 
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C      OFFSPRING  RELATIONSHIPS  ARF  KNOWN*  THUS  PERMITTING  COVA»lA'lCES 

C      AMONG  FAMILIES  OP  THE  CURRFNT  GENERATION  TO  RF  COmPUTFD  «*«-**»*«»* 

C 

K  =  0 

NIOLO=0 

UO  144  I=1»NI 

IF  (I0AM(I)-IUAM(I*1) )  140,144*140 

14ft     K=K*1 

INnWlN(K)=I-NlOLD 

NlOLD=I 

00  141  J=1*NIP 

IF  (ISIRE(I) .NE.ISTORI (J) )  Go  TO  141 

IFAMS(K)«FCOUEI (J) 

GO  TO  142 
141     CONTINUE 

IFAMS(K)=0 
14?     DO  14.1  J=1»NIP 

IF  (TDAM(I) .NE.ISTORT (J) )  GO  TO  143 

IFAMD(K)=FCOUEI (J) 

NEWFAM(K)=KCODEP(T) 

GO  TO  144 
141     CONTINUE 

IFAMD(K)=0 
NEwFAM(K)=FCOUEPU> 

144  CONTINUE 
KCOUNT=K 

c 

C      READ  AND  STORE  RELATED  FAMILIFS  ANp  COVARTANCES  FpOM  PREVIOUS  GFN. 
C 

ND  =  0 

DO  14Q  M=l ,NUMFP1 
I1=M 

READ  (5*147)  L(M) » <IR(I1 »J) *IC<T1*J) *J=1*9) 
IF  (L(M))  150*150*145 

145  NDeND*l 

IF  (L(M)-9)  149,149,14* 
14ft     LVALUE=L(M) 

READ  (5«14R)   UR(IltJ)  ,TCIIliJ)  ,  J=l  0  , LVAl.UF  ) 

149  CONTINUE 

150  WRITE  (6,151)  IGEN,NM,NP 
WRITE  (6,1?4) 

WRITE  (6*152) 

C 

C      COMPUTE  FAMILY  INBREEDING  COEFF I C IfNTS»********«>*********»******** 

C 

DO    lflf,    LS  =  1*KC0UIMT 
K  =  l 

COv/  =  0. 
DO    156    LT=1«NU 

IF  (IR(LTtl)-IFAMS(LS) )  156*153,156 
15^        MATCHS=LT 

LIMIT=L(LT) 

DO  155  INHH=1»LIMIT 

IF  (TFAMU(LS)-IR(LT*1NRR) )  1^5,1S4,155 

154  C0VslC(LT*IN8R) 
TNBRDC(LS)=(CO"/?.) 

IMeSUM=lNpSUM*TNPROC (LS) *TNOWTN (L S) 
INBKDC(LS)=INbRDC(LS)*l000C 

IREL (K)sNEWFAM(LSl 
GO  TO  157 

155  CONTINUE 
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INBRDC(l.S)  =  (AVUROP/?.)«inOOO. 
INHSUM=TNBSUM*lNBPDr (LS) *INDw)N(lS) 
INBPnC(LS)=INBHl)r  (LS)  *10000 
IREL(K)=NE*FAM(LS) 
GO  TO  157 
15ft     COmTINUE 

WRTTF  <6»161>  IFAMS(LS)  ,IFAMD(LS) ,NFWFAM(LS) 
NUmINV=NUMINV*1 
GO  TO  IB? 
C 

C      COMPUTE  COVAfrlANCE  WITHIN  FAMTL  Y*********************************** 
C 

157     HO  160  LU=1»NU 

IF  (IR(LUti)-IFAMO(LS))  160»15B,160 
15R        MATCHD=LU 

IF  (IFAMS(LS)  .NE.IFAMD(LS)  )  G.O  TO  159 
XslNBRDP(LU) 
lCOV(K)=(CUV+X)/2. 
GO  TO  lft2 
150        X«lNRROP(LT)*lNBRnP(L'J) 
IrOV(K)=(CUV*X/2.)/?. 
GO  TO  16? 
16n     CONTINUE 

WRTTE  (6»161)  IFAMS(LS)  ,IFflMH(L«;)  ,NFWFAM(LS) 
NUmINV=NUMTNV*1 
GO  TO  1R5 
C 
C      SEARCH  FOR  RtLATtO  FAMILIES  AND  COVARIANCFS, PUNCH  AND  PRINT  OUTPUT 


c 


16?     LIMITS=L(MATCHS) 
LIMITD=L(MATCHD) 
DO  lfll  N=1»KC0UNT 
COV=0. 
DO  169  LlStHS=ltLIMITS 

IF  (IFAMS(N)-IR(MATCHS,LISTRS) )  166«163»166 
16?  TF  (NEWPAM(N)-NEWFAM(LS) )  I64»16ft«l64 

16*  CV=IC(MATCHS,LISTRS) 

NFND1 0=NEWFAM(N) /\0 
MFLD10=NEWFAM(LS)/10 
IF  (NFNU10.NE.NI-LU10)  GO  TO  165 
X=INBRDH(MATCHS) 
COV«COV*X 

GO  TO  lb* 

165  COV=COV*CV 

166  IF  (IFAMU(N)-IR(MATCHS,LISTR«;)  )  169»l67»l69 

167  IF  (NEWF Am(N)-NFWFAM(LS) )  16«.169»l6H 
16P            CV=IC(MATCHS,LTSTRS) 

COV=COV*CV 
16Q        CONTINUE 

DO  175  LISIRO=l.LTMlTn 

IF  (IFAMS(N)-IR(MATCHD,LI<;TRnn  17J>tl70»l72 
170  IF  (NEWKAM(N)-NEWFAM(LS) )  I71tl7?tl71 

17)  CV=IC(MATCHUtLISTRD) 

COV=COV*CV 
17?  IF  (IFAMD(N)-IR(MATCHD,LI«;TRn)  )  175»173,l75 

17?  IF  (NEWFAM(N)-NEWFAM(LS)  )  l74»17=;»l74 

174  CV=IC(MATCHD»LTSTRD) 
COV  =  COV*Cv/ 

175  CONTINUE 

IF  (COV)  1?6»181i176 
17ft        IF  (COV.LE.COVUEF)  GO  TO  1B0 
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K  =  K*1 

IF  (K-NUMRtL)  179, 179*177 
177        WRITE  (6,176)  NUMREL ,NEWFAM (N) 
KsK-1 
GO  TO  1R2 

179  ICOV(K)=COV/4. 
IREL (K)bNEWFAM(N) 
GO  TO  lfll 

180  COVACC«COVACC*COV 
ACCNUM  =  ACCNUM.*1. 

181  CONTINUE 

187     IF  (K-9)  183,183*184 

183  PUNCH  134*  K» (IREL(M) fICOV(M) iMsl,K) 
WRTTE  (6,126)  Kt  (IREL(M)  ,IC0V(M)  ,Msl,K) 
KArCUM=KACCUM*K 

GO  TO  185 

184  PUMCH  134,  K,  (IPEL(M)  ,ICOV(M)  ,M-1,Q) 
PUNCH  135»  (IKEL(M)  ,ICO\/(M)  ,M=l0»^) 
WRITE  (6,126)  K,  (IREL  (M)  ,ICOV(M)  ,M=1,9) 
WRITE  (6*125)  (IREL(M) ,TCOV(M) tMxio,K) 
KArCUM=KACCUM*K 

185  CONTINUE 
PUNCH  119 

PUNCH  139,  (INBHUC(J) »J=1,NFC) 
C 
C      CONVERT  COVARTANCES  OF  INDIVIDUALS  WITHTN  FAMILIES  TO  INBR.  CnEFF. 


C 


DO  18ft  JPxl,NFC 

C0FFF(JP)=TNBHDC(JP) 
COFFF( jP)=(COhFF(JP)/l0000.)-l. 
18ft  CONTINUE 

WRITE  (6,187)  (CUEFF(J) ,J=l,NFC) 

UsNM-NUMINV 

IF  (O.LT.l.)  0=1. 

TOTINRslNRSUM 

AVGINR=(TOTINR/D)/100. 

TOTINRcTOTlNB/lOyOO. 

01=NFr 

AVNOORsACCNUM/Dl 

D2sKArCUM 

AVNORl  a(D2*ACCNUM) /Ol 

IF  (ArCNUM.GT.O.)  ^0  TO  188 

ACCNUmsI, 
18ft  AVGDRPs(C0VACC/4.)/ACCNUM 

AVGDRp=AVG0RP/ 10000. 

WRITE  (6,1R9) 

WRITE  (6,190)  IGtN,NM,MALFS,NFC,NUMlNV,AVN0RL,AVNODR,flV0ROP,AvGnRP 
IiAVGImB 
191  STOP 

101  FORMAT  (5I3»F4,0»F4.4) 

104     FOPMAT  (4X14,2X14,2X14) 

10ft  FORMAT  (Aft) 

107  FORMAT  (  55H  PAHAM^TEW  CARn  NO.  2  OMITTEO  OR  INCORRECT  VALUE  FOR  N 

1M) 
111     FORMAT  (/5H  SlRE»I6»2feH  MATEO  Tn  MORF  THAN  9  OaMS/) 
117  FORMAT  (IS,  27H  PROGENY  COOES  DERlvEO  FR0M,j4.  21M  FAMILIFS  OF  GFM 

1.  NO. ,13) 
lift  FORMAT  (10(214) ) 

119  FORMAT  (1H  ) 

120  FORMAT  (1X15,  27H  PROGENY  CODES  DEPlVFO  FR0m,I4,  ?1h  FAMILIES  OF  G 
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lEN.    NO. ,13) 
121     FORMAT     (132H  LODE  CODE  CODE  CODE  COO 

IE     rOOE     CODE     CODE     CODE     CODE     COO 

2E     COD  E»/) 
12?  FORMAT  (12(1X215) ) 
124  FORMAT  (  B7H  CURkEMT  GENERATION  RELATFO  FAMILIES  AND  COVARIANCFS  W 

1ITHIN  AND  RETWEEN  FAMILIES  FOLLOW,/) 
12*  FORMAT  (1«X9(I5»X5) ) 
12ft  FORMAT  (/ilSt l3X9(I5tI5) ) 

134  FORMAT  (I?,9(T3.Ib) ) 

135  FORMAT  (2X9(l3»Ib)  ) 
137  FORMAT  (15*27X14) 
13ft  FORMAT  (10(214) ) 
130  FORMAT  (1615) 

147  FORMAT  (I2f9 ( X3i 15) ) 
14*  FORMAT  (2X,9(I3,Ib) ) 
151  FORMAT  (/♦  21H  CURRENT  GENERATION  *,T5,/,  2lH  No.  MASTER  CARDS    = 

1 f I 5 ♦ / ,  18H  PUNCH  OUTPUT  PROM, 15,  55H  FAMILIES  IN  PREVIOUS  GtN,  HAS 

2  BEpN  READ  IN  ANU  STORED.,/) 
15?  FORMAT  (/,108H  NO.  REL.  FAM.      FAM,  COV.  FAN.  COV.  TAM.  COV.  FAM 

1.  COV.  FAM.  COV.  FAM.  COV.  FAM.  COV.  FAM.  COV.  FAM.  COV.*/) 
161     FORMAT  (  64H  INVAL.  PFDTGREF-PHF  V I  OIlS  GFN.  PROGFNY  C 

1  ODF  OF  SIRE  =»I5»   8H  OF  DAM=,l5»  35H  CURRENT  REN.  PROGF.NY  COOF 

2  Or  IN0.=*I4,/) 

17ft        FORMAT  (  4bH  NIJMBFR  OF  RFLATFD  FORTES  EXCEEDS, T 

1       a,  38H  PROGENY  CODE  OF  NEXT  RELATED  FAMILY  =»I5,/) 

187  FORMAT  (/,  77H  THE  FOLLOWING  ARE  FAMILY  INBREEDING  COEFFICIENTS  FO 
1R  THE  CURRENT  GENERAT ION. ,/ ( 16F6.4) ) 

189  FORMAT  (//♦  57H1 SUMMARY  FRF SUMMARY  FILE 

1 »//) 

190  FORMAT  (  18H  GENERATION  NO.   =,I7»//,  18H  NO.  INDIVIDUALS  =.l7»//, 

1  18H  MO,  OF  SIRES     =»l7,//»  1RH  NO,  OF  OAMS      =*l7»//»  18H  NO. 

2  INVAL.  PFD.  a, 17,//*  4QH  AVERAGE  NUMp.FR  OF  RELATFD  FAMILIES  IN  TH 
3IS  GEN.=»F6.1*//»  49H  NO.  OF  DROPPED  COVARI ANCES-.CURPFNT  GENFRATT 
40N  s,F6.1,//»  49H  AVERAGE  DROPPED  COVARI ANCE-PKEVTOuS  GENERATION  = 
5tF8,4,//,  49H  AVERAGE  DROPPED  COVARI ANCE— CURRENT  GENFRATION  =,F8. 
64,//,  49H  AVERAGE  PERCENT  INRRFEDlMG — CURRENT  GfnfRATTON  =,FB.4,/) 

END 
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